Spin wave modes in magnetic nanodisks under in-plane magnetic field 
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The size dependence of spin wave modes in a circular Permalloy (Py) nanodisk under an in-plane 
magnetic field is systematically studied by using micromagnetics simulations. We show that as the 
disk diameter is increased, the resonance frequency of the backward mode deceases while that of 
the uniform mode increases. The avoided crossing of resonance frequencies of the uniform mode 
and the backward mode appears in the plot of the size dependence of resonance frequencies and the 
backward mode turns into the so-called "edge mode" for large nanodisks. 
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The dynamics of magnetizations confined in a mag- 
netic nanostructure is of special interest due to emerging 
applications in spintronics devices, such as future record- 
ing head sensors, magnetic random-access memories and 
spin torque oscillators^i Ferromagnetic resonance (FMR) 
measurement is a powerful tool to investigate the high- 
frequency response of magnetizations. When magnetiza- 
tion is confined in a nanoscale structure, the FMR spec- 
tra of the magnetic nanostructures show multi-peaks, i.e., 
several spin wave modes are excitedj^i^ 

The spin wave modes of in-plane magnetized film are 
classified into two modes depending on the relative di- 
rection of the magnetization vector at equilibrium, Mq, 
and the wavevector of the spin wave, qA^ One is the Da- 
mon and Eshbach (DE) mode where q is perpendicular 
to MqA- The other is the backward (BA) mode where q is 
parallel to iWo.— These two modes have been observed in 
magnetic wires using Brillouin light scattering(BLS)i^ 
For magnetic nanostructures, another spin wave mode 
called the "edge mode", where the spin wave is local- 
ized around edges, appears. The edge mode has been 
observed in thin films f^ii^ wires, circular disksiSrJ^ and 
elliptic disks. ^''"^^ It is important to study the relation 
among these spin wave modes in magnetic nanostructures 
from both scientific and practical points of view. 

In this paper, we systematically investigated the size 
dependence of spin wave modes in the magnetic nanodisk 
shown in Fig.[l](a) by using micromagnetics simulations. 
Until now it has been widely accepted that the origin of 
the edge mode is the pinning of the spin wave due to 
the strong demagnetization field around the edges and 
that the edge mode has no relation with the BA or DE 
modes.--'- However, we show that the edge mode can be 
regarded as the BA mode with two nodes. 

In order to solve the Landau-Lifshitz-Gilbert (LLG) 
equation for a circular Py disk with thickness h = 10 
nm under an in-plane magnetic field we employed nmag 
codei^ii^ which is a micromagnetics simulation pack- 
age based on the finite-element and boundary-element 
methods. We assumed a saturation magnetization of 
Ms = 0.8 X 10^ A/m, an exchange stiffness is A = 
1.3 X 10~^^ J/m and a Gilbert dumping constant of 
a = 0.01. The finite element mesh was generated us- 



ing Gmshj2£ The size of the mesh was about 5 nm and we 
confirmed that this mesh was sufficiently fine to describe 
the spin wave dynamics of the magnetic nanodisk we con- 
sidered. The static magnetic field, ffcxt = H^^t^y, was 
applied in the y-direction and the oscillating magnetic 
field, ^rf = hj-f sm{2Tr ft)ex, of frequency / was applied 
in the x-direction as shown in Fig. [1] (a). We assumed 
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FIG. 1. (Color online) (a) Schematic illustration of a circular 
Py nanodisk under static iJcxt and oscillating hrt magnetic 
fields, (b) Real and imaginary parts of the calculated com- 
plex susceptibilities are plotted by the dotted and solid lines, 
respectively, for d=50, 100 and 200 nm. 
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FIG. 2. (Color online) (a) Directions of the static and oscillating magnetic fields are sliown. (b) Snapshots of the distribution 
of the deviation, m{r,t), in a nanodisk with d = 30 nm (i), 50 nm (ii), 70 nm (iii), 90 nm (iv), 110 nm (v), 130 nm (vi), 150 
nm (vii), 170 nm (viii), 190 nm (ix), 210 nm (x), 230 nm (xi) and 250 nm (xii). The pink arrows indicate the directions of the 
deviation, m{r,t), at the anti-nodes. The labels BA-n and DE-n stand for the BA and DE modes with n nodal surfaces. 



Hcxt = 1-0 kOe and h^^ = 1 Oe. We obtained a single do- 
main state as an equilibrium magnetization configuration 
for nanodisks with d < 250 nm. Starting with the equilib- 
rium magnetization configuration we calculated the time 
evolution of the magnetization M{r,t) for 20 ns under 
the in-plane magnetic field of iJoxt + ^rf(i)- 

We introduce the deviation of the magnetization de- 
fined as m{r,t) = M{r,t) — Mo{r), where Mo{r) is 
the equilibrium magnetization. The deviation m{r,t) 
represents the amplitude of a spin wave excited by an 
oscillating magnetic field. The real and imaginary part 
of the complex susceptibility x was obtained by fitting 
the average of the a;-component of the deviation, nix , to 

nix = Re(x)/irf sin(27r/t) - Im(x)/irt cos{2n ft), (1) 

after the oscillation was stabilized. In Fig. [T] (b) we plot 
the real and imaginary part of x by the solid and dot- 
ted lines, respectively for d—50, 100 and 200 nm. The 
resonance frequencies were 6.8 GHz for d = 50 nm, 7.2 
and 10.7 GHz for d = 100 nm and 6.5, 9.3, 11.0 and 
13.8 GHz for d = 200 nm. For d = 200 nm, the mag- 
nitudes of X at 11.0 and 13.8 GHz were so small that 
we used log-scale plots to identify them. The calculated 
resonance frequencies were in good agreement with those 
observed in TR-MOKE measurement of circular Py nan- 
odisks with h =10 nm by Shaw et al. They observed 
about 7 GHz for d = 50 nm, 7 and 10 GHz, for d = 100 
nm, and 7 and 9 GHz for d — 200 nm under the in-plane 
magnetic field of 1 kOe^'^. 

Figure[2] (b) shows snapshots of the deviation, Tn(r, t), 
at the resonance frequencies for d = 30, 50, • • • , 250 nm. 
The directions of the static and oscillating magnetic fields 
in this plot are shown in Fig. [5] (a). The labels BA-n 
and DE-n stand for the BA and DE modes with n nodal 
surfaces as shown in Fig.[5](b). The pink arrows indicate 



the directions of the deviation, Tn(r, t), at the anti-nodes. 
We could not identify the third lowest mode for d = 250 
nm since its resonance frequency was very close to that 
of the center mode. For d = 30 nm, the deviation, m, 
was uniformly distributed throughout the disk, i.e., the 
uniform mode was excited. As the disk diameter, d, 
is increased, other resonant spin wave modes classified 
as BA-2, BA-4 and DE-2 modes appear. It should be 
noted that for a large nanodisk the lowest mode is not the 
uniform mode but the edge mode, and the second lowest 
mode is not the BA-2 mode but the center mode. Tracing 
the spin wave modes of the lowest-resonance frequencies, 
i.e., the bottom line of Fig. [2] (b), one might come up 
with the idea that the origin of the edge mode is the 
uniform mode and it has no relation with either the BA 
or DE modes. However, our systematic investigation of 
the resonant spin wave modes revealed that the above 
interpretation is not valid. 

In Fig. [3] we plot the size dependence of the calculated 
resonance frequencies by circles with radii proportional 
to the imaginary part of the susceptibility, Im(x). The 
blue dotted lines are visual guides connecting resonance 
frequencies in the same branches. The purple solid line 
represents the resonance frequencies given by the Kittel 
equationSi: 

/k = ^fHUH + {N,-Ny)fM). (2) 

where fn = 7i/oxt/27r, Jm = -iMs/^ir, 7 is the gyromag- 
netic ratio and A^'s represent the demagnetization factors 
estimated by means of the micromagnetics simulations. 
For small nanodisks with d < 60 nm, we obtained a single 
resonance peak at frequencies close to those given by the 
Kittel equation of Eq. As the disk diameter, d, is in- 
creased, the frequency of the uniform mode increases due 
to the increase of the demagnetization field, while those 
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FIG. 3. (Color online) Calculated resonance frequencies, /res, 
are plotted by circles with radii proportional to imaginary 
part of the susceptibility, Im(x), against the disk diameter, d. 
The blue dotted lines are visual guides connecting resonance 
frequencies in the same branches. The purple solid line repre- 
sents the resonance frequencies given by the Kittel equation 
of Eq. ((2)1 . The labels of the excited spin wave modes are the 
same as those in Fig. [5] (b). 



of the BA-2, BA-4 and DE-2 modes decrease. At around 
d ~ 150 nm where the frequency of the BA-2 mode comes 
close to that of the uniform mode, an avoided crossing of 
the BA-2 mode and uniform mode appears and the BA-2 
mode becomes the edge mode for large disks as indicated 
by the black dashed line. Therefore, the edge mode can 
be regarded as the BA-2 mode with two nodes. The uni- 
form mode becomes the center mode, the frequency of 
which is close to that given by the Kittel equation. 

In summary, we investigated the size dependence of 
spin wave modes in a circular Py nanodisk with d < 
250 nm under an in-plane magnetic field using micro- 
magnetics simulations. We showed that the excited spin 
wave modes are classified into uniform mode, edge mode, 
center mode, BA-2 mode, BA-4 mode, and DE-2 mode 
depending on the size of the nanodisk. For large disks, 
we found an avoided crossing of the BA-2 and uniform 
modes, at which BA-2 mode (uniform mode) becomes 
the edge mode (center mode). 
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